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I. INTRODUCTION 

It is commonly agreed that at high spacetime curvatures, when the quantum effects of the grav- 
itational field become significant, General Relativity needs to be replaced by a quantum theory of 
gravity. In spite of progress in several directions, finding this new theory has proven a challeng- 
ing problem for several decades. Current research in the field is paying substantial attention to 
the numerous indications that gravity may only be emergent, meaning that it is a collective, or 
thermodynamical, description of microscopic physics in which we do not encounter geometric or 
gravitational degree of freedom. An analogy to illustrate this point of view is fluid dynamics and 
the transition from from thermodynamics to the kinetic theory. What we currently know is the 
low energy theory, the analogue of fluid dynamics. We are looking for the microscopic theory, the 
analogue of the quantum molecular dynamics. Just as there are no waves in the molecular theory, 
we may not find geometric degree of freedom in the fundamental theory. Not surprisingly, this 
significant shift in perspective opens up new routes that may take us out of the old problems. 

The emergent viewpoint amounts to treating quantum gravity as a problem in statistical physics. 
A powerful set of methods in statistical physics involve the use of lattice-based models, such as 
the Ising model for ferromagnetism, the Hubbard model for the conductor /insulator transition, 
etc. Such methods are starting to be introduced in quantum gravity. Examples are G. Volovik's 
work on emergent Lorentz invariance at the Fermi point [TJ, X.-G. Wen's work on emergent matter 
and gravitons from a bosonic spin system [2], the emergence of a Lorentzian metric and aspects 
of gravitation such as Hawking radiation in analog models of gravity [3J, as well as long-standing 
approaches such as matrix models [3], and more radical formulations of geometry in terms of 
information [5j [6] . 

It is natural for the lattice in the lattice system to play the role of (a primitive form of) 
geometry. Now, General Relativity is a background independent theory, by which we mean that 
the geometry of spacetime is fully dynamical. By analogy, we expect that the use of a fixed lattice is 
inappropriate and one instead needs models on a dynamical lattice. While for this reason desirable, 
dynamical lattices raise difficult technical problems that have not been previously addressed in the 
field of statistical physics. The present article is concerned with exactly this problem and presents 
a method that deals with dynamical lattices in certain situations. 

In previous work, we introduced Quantum Graphity, a background independent model of space- 
time in which time is an external parameter and space is described by a relational theory based 
on graphs [13]. The idea is to represent locality by the adjacencies of a dynamical graph on which 
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the diffeomorphism group is replaced, in the high energy phase, by the symmetric group on the 
complete graph which breaks down to its subgraphs at lower energies. These kind of models are 
sometimes called event-symmetric^^) . The principle of event symmetry refers to the replacement, 
at high energies, of the diffeomorphism group with the group of permutations of events in space- 
time. In the context of Quantum Graphity, however, the event symmetry is only spatial, in the 
sense that at high energy the graph is complete and every vertex of the graph is at distance one 
from each other. If the dynamics is such that the system settles into a minimum energy subgraph 
that exhibits geometric symmetries, for instance, a discrete version of flat space in low dimension, 
we say that geometry emerges in that phase. In |13| was shown that desired symmetric lattices, 
i.e., discrete 2d FRW, are stable local minima under certain choices of parameters in the hamilto- 
nian (when the effect of the matter on the lattice is neglected). Following this work, in |16j . we 
used the same concept of locality in terms of a dynamical lattice, but with a new type of matter 
that interacts non-linearly with the geometry, a precursor of gravity, and initiated a study of the 
quantum properties of that system. 

In the present article, we return to the technical issues of spin systems on dynamical lattices 
and we introduce a method to deal with a theory of dynamical graphs on N vertices. Such graphs 
are subgraphs of IK at, the complete graph on N vertices. We show how, by transforming %n to 
its line graph, the theory can be approximated by an Ising model on the line graph of a complete 
graph. We then use this to study the low energy properties of the Quantum Graphity model in 
|13j . Using mean field theory we calculate the average valence of the graph at low temperature 
and we evaluate the first corrections due to the presence of 3- loops. 

The paper is organized as follows. In section II, we review the Quantum Graphity model [13] 
with no matter. In section III, we define the line graph derived from a generic graph and summarize 
its properties. In section IV, we show how, in a certain reasonable approximation, the hamiltonian 
of Quantum Graphity can be recast as an Ising hamiltonian on the line graph of a complete graph. 
In section V, we calculate the corrections due to loops at low temperature and describe, in this 
framework, the behavior of the correlation function in mean field theory. Conclusions follow in 
Section VI. 

II. QUANTUM GRAPHITY 

Let us briefly introduce the Quantum Graphity model [13J. As the name graphity implies, 
Quantum Graphity is a model for a quantum theory of gravity in which the fundamental microstates 
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are dynamical graphs, postulated to describe relational physics at Planckian energies. There is no 
notion of geometry or quantum geometry at high energy, instead, geometry emerges as the system 
cools down and away from the Planckian regime. The microstates live in a Hilbert space on the 

complete graph X/y of N vertices, given by 

MA^i, n 

U 3 

where en represents the edge of the graph Kn between the i and j vertices, while and IK 1 / are 
the Hilbert spaces associated with edges and vertices respectively. In particular, the Hilbert space 
associated with an edge between vertex i and j is the two-level state space: 

^.=span{|0),|l)}. (1) 

The two states 1,0 in ([!]) are interpreted as the edge being on or off respectively. This choice 
means that basis states in the Hilbert space of the edges represent subgraphs $x N of the complete 
graph %n. A generic state in the Hilbert space of the edges is a superposition of such subgraphs: 

1^) = ^2c t 9K N ;t- 

In the full model of [T3], extra degrees of freedom are assigned to the on states: 

5C£ = span{|0)y, \h)ij, \h)ij}- (2) 

In |13| . and in the present work, there are no degrees of freedom associated with the vertices and 
hence we ignore JC". 

Let us focus now on the state space 0. On the Fock space of the edges we can introduce the 
ladder operators a and a), with the usual action: 

a\j\l)ij = ciij\0}ij = , a\l)ij = \0)ij. 

Dynamics in Quantum Graphity is given by a hamiltonian acting on the graph states of the form[13j: 

H = Hy + Hl + interaction terms, (3) 

where Hy keeps track of how many on edges are attached to a single vertex, and Hl counts closed 
paths in the graph. The interaction term will not be used in the following, but in a generic model 
these terms produce Alexander moves on the graph. 

In more detail, the term Hy is usually chosen to be of the form: 

j 
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where the indices i,j = 1, ./V enumerate vertices, iVy = aLay is the usual number operator, and 
gv and p are free couplings that we assume to be positive. The purpose of this term is to ensure 
that at low energies the system has a (low temperature) phase in which the average vertex valence 
(i.e. on edges attached to a vertex) is vq. Later on in the paper we will show that, at least in the 
mean field theory approximation, this is indeed the case. The term Hl is given by 

i L 

where g^ and r are couplings assumed to be positive. The operator P(i,L) counts the number 
of non- retracing paths of length L based at the vertex i. This operator is related to the trace of 
the adjacency matrix in the original model. We will build this operator in another way in the 
following. For r < 1, so that higher length loops contribute less than short length loops, this term 
is semi-local. The L\ comes from the expansion of the exponential of the loop-path operator [13J. 
It was shown in [13] that r determines the length on which loop size is peaked at low energies. 

In what follows we introduce a new method to analyze dynamical lattices, by transforming %n 
to its line graph which we define in the next section. The on/ 'off edges of %n will become Ising 
spins on the fixed line graph, so that standard methods of statistical physics can be used. 

III. GRAPHS AND LINE GRAPHS 

We start by defining line graphs. Let G = (V, E) denote a graph with vertex set V = {v\, i»2, •••} 
and edge set E = {ei,e2, ...}. The line graph £>(G) = (V,E) is the graph of the adjacencies of 
G, containing information on the connectivity of the original graph. Each vertex v G V{L{G)) 
corresponds to an edge e G E(G). Two vertices v\ and V2 in V(£(G)) are adjacent if and only if 
the edges in G corresponding to v\ and V2 share a vertex. The correspondence between G and £(G) 
is not one to one. From a given graph G we can construct only one £(G) but it is not true that 
any graph is a line graph. In fact, according to the Beineke classification, there are 9 non- minimal 
graphs that are not line graphs of another graph and each graph containing them is not a line 
graph[T7]. The simplest example of a line graph is depicted in Fig. [T] 

Given a graph G, we can construct its line graph using the following procedure: 

1. Enumerate the vertices of G. 

2. Enumerate the edges of G with a fixed prescription (see example below) and put a blob on 
them. 
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FIG. 1. The simplest example of element of the Beineke classification of a minimal graph that is not a line 
graph of any other one. This means that in general there is not a one-one correspondence between a graphs 
and line- duals. 

3. If two edges share a vertex, draw a bold line between them. 

4. Remove G and its enumeration. 

What is left is the line graph of G where the blobs represent its vertices. 
Let us now introduce some useful quantities: 

Definition. (Kirchhoff matrix) Let G be a generic graph, V = {v\, V2, • • • , v n } be the set of vertices 
of G and E = e2, • • • , e p } be the set of edges of G. Let us define the matrix P of size n x p 
with entries Pip, where i is an integer between 1 and n on the set of vertices and f3 is an integer 
between 1 and p on the set of edges, such that: 

{1 if the edqe B has an vertex on the vertex i, 
(6) 
otherwise . 

The Kirchhoff matrix K is the p x p matrix built from P, such that: 

K = P l P, (7) 

pt re p resen ti n g the transpose of P. 

A well-known theorem now gives the incidence matrix of the line graph L(G): 

Theorem. Let G be a graph with p edges and n vertices and let £(G) be its line graph. Then the 
matrix: 

J = K — 2 I, (8) 

where I is the p x p identity matrix, is the incidence matrix of L{G). 

In the next section we will show how the graphity hamiltonian can be recast on the line graph 
using (|6j) and 
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IV. THE LINE GRAPH REPRESENTATION 

Since in the hamiltonian ^ we neglect the terms in which vertices are interacting because we 
assume there are no degrees of freedom on them, one could expect that it can be rewritten only in 
terms of the connectivity of the graph. To carry out such a reformulation, let us expand the first 
term in ^ for small values of the parameter p: 

Hv = gv^i+pgv 

= g v (1 + vlp) 1 +P9V E fiiifyk ~ ^9vpv E #v + °(P 2 )- ( 9 ) 

i ijk ij 

As we will see later, such an expansion does not modify the properties of the model at low tem- 
perature. The first term in ^ is an energy shift and, for what is to come, can be neglected. We 
should now be able to recognize some particular terms in the expansion ([9]). The third term is 
proportional to the operator Y^ij ^ij- It is the sum over all the edges of the graph, zero or not, of 
the number operator. We will change the notation to 

E^^ 2 E% 

ij P 

where /3, as in the previous section, runs from 1 to N{N — l)/2 and labels the edges of %m or, 
equivalently in what follows, the vertices of its line graph. 

To rewrite the second term in ([9]), we need the matrix Pip of ([6]) in this context. This matrix 
maps the graph to its line graph, as we will see. We first fix a prescription to label edges. Let %n be 
the complete graph of N vertices. Let J be any enumeration of V(%n), i E J = 1, N. We identify 
edges by their endpoint vertices (i, j), with i, j £ J. A labeling Sp, j3 € £ = {1, • • • , ./V(iV — l)/2} 
is an enumeration of E(L(%n)), according to the following prescription: 

Si,- - - , SV-i label the edges connecting the vertices {(1, 2), • • • , (1, N)}; 
Sn, ■ ■ ■ ,S 2 w-i) label the edges connecting the vertices {(2,3), • • • , (2, N)}; 



Sjv(7v_i)/2 labels the edge connecting the vertices (N — 1, N). (10) 

Using this prescription it is easy to see that the matrix introduced in Q, for the complete 
graph %n, has the simple (recursive) form: 

P»=\^ 6 |, (11) 
1 tJV— 1 pN-l ' 
L b'd r a'a.' 



s 




FIG. 2. Three examples of complete graphs labeled according to the prescription (10 1 



where Vn-i is a row vector of length N — 1, I^ -1 is the identity matrix of size (N — 1) x (N — 1) 
and represents a null row vector of length N(N — l)/2 — (N — 1). The indices {a', b', c'} and a' 
run from 1 to iV — 1 and 1 to (N — 1)(N — 2)/2 respectively. As an example, for the graphs of Fig. 
[2] the P matrices are: 



fl 1 o\ 
1 1 

V° 1 l J 



(12) 



( 1 


1 


1 










1 








1 


1 








1 





1 





1 


v° 





1 





1 


1 J 



(13) 



/l 


1 


1 


1 

















o\ 


1 











1 


1 


1 














1 








1 








1 


1 











1 








1 





1 





1 


1° 








1 








1 





1 


1 ) 



for (a), (b) and (c) respectively. 

It is easy to see that two edges a and f3 have a common vertex if and only if we have: 



(15) 
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(16) 



where * is the transposition operation. By construction, c can take the following values only: 

2 if a = 0, 

1 if a /3 and a and /3 have a common vertex, 
if a 7^ (3 and a and /3 do not have a common vertex. 
In particular, for TV = 4, K is given by: 

/2 1 1 1 1 (A 
12 10 1 
112 11 
1 2 1 1 
10 12 1 
\ 1 1 1 1 2 J 

Using now the matrix P^p just introduced, we want to construct generic n-string matrices as 
composition of n edges of the graph, thus a path on the graph. The n-string matrices will be 
needed both for the non-retracing loop term and the 2-edge interaction term in equation ([9]). The 
quantity: 



K 4 



(17) 



(18) 



is the definition of the Kirchhoff matrix of equation (|8|) if we sum over the index i. From K % aB we 
can construct strings of P's of the form 



Q 



K n K 12 ■ ■ ■ K ln 



that we call string matrices of nth order. These string matrices represent paths through vertices 
i\ - ■ - in and they are zero unless the edges corresponding to a\ ■ ■ ■ a n+ \ are in the correct order, 
that means, they represent an actual path on the graph. For instance, the number of paths of 
length 2 on the complete graph is given by 



#2-strings= ^ K *fl, 



or, equivalently, we can use (|6j) and rewrite @ as 



#2-strings = Y.^= ^ ~ 2 ^ 



(19) 



(20) 



The subtraction of twice the identity in ( 20 ) is the same as the subtraction of the self-energy of 



each edge. We now clearly see that this matrix is precisely the incidence matrix of the line graph 
ofX N introduced in (§), with a,p € V{L(X N )). 
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So far we have dealt with the complete graph only. We wish to extend this formalism to a 
dynamical graph. In order to do that we return to the Hilbert space formulation of the graph with 
on/off edges. Recall that any graph on N vertices is a subgraph of the complete graph %n, with 
some edges off. Thus, since we can always map a graph on a complete graph, we can count paths 



on any graph by modifying (20) so that it counts paths of only on edges on the corresponding 
complete graph. To do so, we introduce in the sum the number operators Ng in the following way: 

#2-strings = ^ (K a3 - 2 I a p) N a Np 

This term does not contribute if any of the two edges a, /3 is off. It is easy to see that this term of 
the hamiltonian is an Ising interaction. The important difference between these two hamiltonians 
is that in our case the spin system is on the line graph a complete graph 3Cjv- 

By extension of the above, we are now able to construct a generic path operator out of K l a oS. 
We define 



Q ai—a n 

= E E N ai ..-N an , (22) 



Q oi-o„ 



where the set Q is 



il ^ ■ ■ ■ 7^ i n -% G J for non-retracing paths, 
2 = { (23) 
ii, . . . , i n —\ £ ? for retracing paths. 

It is easy to see that it counts the number of paths of length n in the graph, that is why we call 

the Q's string matrices. Note that Qai'-Sn 1 can take values and 1 only because it is a product of 

0's and l's. This string matrix is not the matrix multiplication of the Kirchhoff matrices: it only 

reduces to matrix multiplication for retracing paths where we sum over all possible vertices. 



In the following, we will denote the two sets in (23) as Q r and Q nr for the retracing and non 
retracing cases respectively; moreover, we may explicitly show the indices on which we are doing 
the sum as Q r / nr {i b ^). In order to count loops, we just need to impose a± = a n : 

p.p. p.p. p.p. 

K h K i2 K is 

Jl aitt2 Jl «203 il Q!3Q!l 

Thus we have discovered that, when there are no degrees of freedom on the vertices of the graph 
and we neglect the interaction terms, we can recast the Quantum Graphity hamiltonian on the 
line graph L(%n) representation in the weak coupling regime at finite N. 
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We end this section with two properties of the nth-order string matrices. Let us define: 

fjr/nr _ ST^ rfii—in-i 

Qr / nr 



(24) 



The following properties of the sum of these string matrices on complete graphs will be required 
next: 

Property 1: Let S = %n- Then, for a loop of n edges: 



E QZ 



■a L ai 



N(N — 1)---(N-L) 



(25) 



and 



E Wax- 



■a L ai 



N 



(26) 



Proof. These two facts follow trivially if we note that the equations ( 25 ) and ( 26 ) count the number 



of retracing and non-retracing paths of length L on the complete graph respectively. 



Property 2: Let S = %n- Then, for a loop of n edges, and for L > 4, we have: 



E^M^-^ ®Z ckc* = (N - 3) • • • (N - 3 - (L - 4))K aia2 , 



(27) 



while, for L = 3: 



Er\nr 



a 2 ■ 



(28) 



if ai ^ a 2 7^ Qi3 ■ • • j=- a n . 



Proof. Note that Yla s ^a 4 =/=-^a n Qcn a n a Y i s the number of non-retracing loops of length L 

on the complete graph %n which pass by the edges a\ and a 2 - Now, it is easy to see that if the 
edges a\ and a 2 do not share a link this quantity is zero. Also note that, by the symmetry of the 
complete graph, the number of non-retracing loops based on two neighboring edges must be the 
same for each pair of edg sharing a node. Since the matrix K a . Q . takes values 1 or 

depending on whether the edges aj 1 ,aj 2 are neighbors or not, J2a 3 ^a 4 ^--^a n Qai uiai mus t De 

proportional to the matrix K aia2 . In order to evaluate the proportionality constant, let us note 
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that each loop is weighed by a factor of 1 because n-string matrices take values 1 or only. The 
combinatorial quantity (N — 3) • • • (N — 3 — (L — 4)) is then the number of non-retracing loops of 
length L passing from two consecutive fixed edges on the complete graph of N vertices, as can be 



easily checked. The special case (29) follows from the fact that if we fix two edges there is only 
one edge which closes the 3-loop. 

Note that, for N L, we have: 

E QZ a^" E %i « Lai =N L . (29) 

We can now collect the results of this section to write the hamiltonian ^ as 

H = A ^ JafiN a Np - B^N a -C Yl %8.7 NcNpN^, (30) 

where 

A = p gv, 

B = 2 g v p v , 

C = 9Lj, (31) 

and neglecting higher order loop terms. 



V. MEAN FIELD THEORY APPROXIMATION AND LOW TEMPERATURE 

EXPANSION 

Having rewritten the hamiltonian in an Ising fashion, we now can approach the problem of 
finding a graph observable and its equilibrium distribution using mean field theory. As we will see, 
the natural graph observable to consider is the average valence of the graph. We will assume that 
the system is at equilibrium and we neglect the interaction terms. In this case, it is straightforward 
to use mean field theory analysis |18j . In what follows, we assume units in which the Boltzmann 
constant ks = 1. 

We start by replacing the number operators N a with semi-classical analogs, imposing that their 
expectation value must lie in the interval / = [0, 1]: 

N/3 -> (Np) P = mp, 
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where P is a probability measure of the following form: 

P{mp) = mp5 1>mfj + (1 - mp)8 Q , mij . (32) 

It is easy to see that this probability distribution forces the spin-average to lie in /. Recall that in 
order to obtain the mean field theory distribution we have to extremize the Gibbs functional given 
by 

$[m] = H[m] - ^S[m], (33) 

where (3 = T -1 , H[ m] is the energy and S[m] is the entropy functional. The latter can be written 
as: 

S[ m ] = - J2 niP(m p (i))logP(m p (i)), (34) 

mp={0,l} 

where ni is the degeneracy of the state. 



Case I: Non-degenerate edge states 

In this subsection we focus on the case in which the states on and off are not degenerate, so 
that rii = 1. In the next subsection we will deal with non-degenerate edge states and in particular 
with 3-degenerate on states. 

In the process of extremizing the Gibbs functional we will see how the average valence of the 
graph naturally emerges. We impose: 

d m ^[m] =0. 

Using 

^S[m] = -log(^_) 

and 

d mf} H[m) = A ^ Japma -B-C ^ Q a ^m a m^, 
we find that the distribution for the m a is 
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The solution of this equation gives the equilibrium value of mp once the value of the temperature 
is fixed. 



We now want to write (35) as a function of an average quantity on the graph. Let us first note 



that, in the mean field theory approximation, we have 

^J Q/3 m a = 2 d(T), (36) 

Q 

where d(T) is the mean valence of the graph. The valence d(T) is a good graph observable that 
we can use also as a double check for our procedure since it appears explicitly in the original 
formulation of the hamiltonian and in the low temperature regime must take the value vq. First, 
we note that: 

N on edges Y,i&d(i)/2 _ d(T) 

TTl n 



N(N- 1)/2 N(N- 1)/2 7V-1 

In the first equality, iV on edges is the number of edges of the graph which are in an on state. In the 
second equality, the average valence (the sum over all the local valencies divided by the number of 
vertices) is explicitly written as a temperature dependent quantity. In the third equality we used 
the graph property: 

The most complicated term in the hamiltonian is the 3-loop one. The simplest way to deal with 
it is to use the Ansatz dictated by the mean field theory: 

/, QaP-yrngmp x> £(T)d 2 (T). (37) 

are® 



Let us replace mp with its average value: d(T)/N — 1. Using eq. (25) for non-retracing paths 
and assuming JV > 1 we obtain the dependence on d(T). £(T) is a function of order ~ 1 at low 
temperature, which we assume is dependent on T. Using these approximations we can see that 
d(T) is a natural order parameter for our mean field theory since it is easily recognized as implicitly 
defined in the stable distribution: 



N — 1 

1 + e md(T)A-a%d?(T)-B} ' 



d ( T ) = t a, m 5 • (38) 



Again, in order to double check our procedure, we can ask if such an order parameter behaves 
as expected at low temperature. We must keep in mind that the starting hamiltonian ([3]) was 
constructed in such a way that the average valence at zero temperature was a fixed value of the 
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parameter vq at finite N. We can now use (38) to check if this is the case. To do so, we Taylor 
expand both sides and match the zeroth and first order coefficients on the left and right hand side 
of the equation. That is, we start with the expansion 

d(T) = a + (3T + 0(T 2 ), (39) 

and, for the approximation to be consistent at T = 0, we require analyticity of the order parameter 
(this has to be the case for a finite volume system in ordinary statistical mechanics, which is the 



case for finite N) . We then require that inside the exponential of equation (38) the temperature 
independent terms in the numerator cancel out so that at T = the exponent is well defined. This 
gives the second order equation in a: 

2aA-£(0)Ca 2 = B. 

Now note that, while this equation has two solutions, we need to only look for the one which is 
analytical in the parameters of the model and tends smoothly to the solution a = in the C — > 
limit. This fixes a to the value a, given by 



We can now plug a at T = into ( 38 ) : 

N — 1 



o 



1 + e {2A-i{0)Ca)P ' 



to obtain the value of /3 in (39): 



^ M-tWa ^fir 1 - 1 )- (41) 

It is easy to see that in the limit ./V — > oo we have /3 — > oo, indicating a second-order phase 
transition (a discontinuity in the first derivative of the order parameter). In our case, this happens 
at T = 0, meaning that this transition is not possible because there is no way to cool down the 
system to zero temperature with an external bath. However, we have to remember that we are just 
approximating the real system with a semi-classical analog. We then simply interpret the above 
result as the fact that the system reaches the ground state very quickly when the temperature 
approaches zero. 



It is interesting now to plug in the couplings. Inserting equations (31) into (40), we find that 
at T = 



g L r 3 v 
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FIG. 3. The behavior of log(d(T)) (y-axis) against T (x-axis), for increasing TV. 



Note that, for small values of r, when r 3 <C , we have a = vq, meaning that at low temperature 
the mean degree is the one imposed by the degree term of the hamiltonian, as expected. We can, 
however, see how the 3-loops term contributes to this quantity by a Taylor expansion in r: 



3 pg v 



(43) 



From this expression it is clear that the loop terms are suppressed if gy 3> gi- This is the 
main result derived in this paper using the line graph representation. A plot of the function d(T) 
is shown in Figure 3. 

We now have the tools to calculate the susceptibility function for the theory in the mean field 
theory approximation. Recall that the susceptibility function tells us how the system reacts to 
a variation of the external magnetic field. In our case, the magnetic field is the combination 
2gvpvQ and we note that the parameter vq in the hamiltonian appears only here. We have the 
following analogy: v o represents the external magnetic field, while 2gyp represents the spin-coupling 
combination. 

In order to calculate the susceptibility, we assume that the constant B is site-dependent (i.e. a 
field). Thus, we have 



in \ 1 9 9 jpihi 1 dm P 

{G a p) B = ~ QD _ flD F[h] - 



(44) 



f3dB p dB a 11 (3dB a ~ 
In particular, we are interested in the susceptibility function when T « 0. From the study of it we 
can gain some information about the low energy behavior of the model. We expand the equilibrium 
distribution 



B, 



(45) 
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Using the notation p a p := ^ 7eS Q Q /3 7 ?7i 7 , we obtain 

Bp = Bp-~ = J2 {^f + AJ aP + CpJ\ m a . (46) 
To invert this equation, we approximate p a p by replacing m 7 — > d(T)/(N — 1): 

7 es 7 es 



We can now can use property (29) of the Q matrices to find that the sum of the Q's reduces to the 



incidence matrix of L(Kn). Hence, inverting equation (46), we obtain 

m P=/ J Q-ypB-y, 
7 

where Q^r = (24 + coJ) -1 . and Co is an effective constant in front of the Ising term of the 

p 7 p 

hamiltonian: 

^ d ( T ) f A»S 

co ~ P9V + ffijfj^. (47) 



It is interesting to note that, thanks to property (28), we can sum all the loop terms up to a 
finite number 1 <C L -C in the hamiltonian if we assume the mean field theory approximation. 
Inserting the couplings, we find 

(48) 

It is interesting to note that in the limit r — > 0, or T — > (where d(T) tends to a finite number for 
N 3> 1), this effective constant tends to pgy- We interpret this as the fact that at low temperature 
the loops become less and less important and the model is dominated by the Ising term. In 
particular, since the external "magnetic" field is given by vq and is assumed to be nonzero, it is 
not surprising that at T = the average valence, the equivalent "magnetization" , approaches this 
value. We note that the N — > oo limit does not behave well unless L = 3. Higher loops are highly 
non-local objects. For a given pair of edges, all the L-loops based on these edges span the whole 



graph already at L = 4, while of course this is not the case for 3- loops. As a result, in formula (48) 
there is a factor proportional to N L ~ 3 which is not present at L = 3. 



Case II: Degenerate edge states 



The Quantum Graphity model [13] allows for degenerate on states on the edges or the vertices 
of the graph. Degeneracy of edge states is necessary, for instance, in order to have emergent matter 
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via the string-net condensation mechanism of Levin and Wen. Degeneracy requires modifying our 
calculations above and we will address it in this subsection. 

The first possible generalization of the Quantum Graphity model is to introduce a Hilbert space 
on the edges of the form ([2]): 

JC| = span{|O>0, 111)/?, |1 2 )/?, |1 3 )/?}- 



This changes the degeneracy number in equation (34). With n\ = 3 and no = 1, we obtain 



3 

d mfj S[m] = -\og{-^—)-2. 
' 1 — mp 

The equilibrium distribution solves this equation. If we put Q = exp[f3(d mi3 H[m] — 2)], we have 

me - (I) 1 « , + (^^TW)\ 

W (9 + V3V27 + 4Q3) 5 2333 
obtained from the only real solution of the third order polynomial equation + Quia — 1 = 0. 



Using the same procedure as before, it is easy to see that equation (40) remains unchanged: 
the low energy average valence is the same in both cases. However, the first derivative, that is, 
the coefficient of the T term in the Taylor expansion of the average valence in the temperature, 
changes, so that ^31 > /J^i (with the obvious notation for the two coefficients). This phenomenon 
can be understood using the following argument. At high temperature, the two models behave in 
the same way, forcing the valence to be high. When the temperature drops, d(T) also goes down. 
While in the (1,1) case the phase space of the on edges is the same as that of the off edges, in 
the (3,1) case the system prefers to stay in the on state. Thus, when the temperature decreases 
the system (3,1) is, at first, slowly converging to the ground state, but at T = it is forced to go 
to the ground state. For this reason the function d(T) has a greater derivative near T = in the 
(3,1) case. 



VI. CONCLUSIONS 



In this paper we introduced a technique to map the Quantum Graphity hamiltonian on the 
line graph of a complete graph. This procedure requires the introduction of the Kirchhoff matrix 
of a graph and the n-string matrices related to these. This mapping is general and not specific 
to Quantum Graphity. Using this mapping in a weak coupling approximation of the model, the 
mean field theory approximation and the low temperature expansion, we studied the properties of 
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the model near zero temperature after having identified the average deg ICC BjS £1 order parameter. 
We found that the model is dual to an Ising model with external nonzero magnetic field if we 
neglect the interaction terms due to loops. In particular, we showed that the average valence is 
naturally a good order parameter for the mean field theory approximation and we found, implicitly, 
its average distribution using the mapped hamiltonian and the mean field theory approximation 
for the 3-loop term. We then studied the susceptibility function and showed how the duality with 
the Ising model can help to interpret the results. In particular, the parameter vo plays the role 
of the external magnetic field. Since vq is assumed to be never zero, the model has no phase 
transition and at T = the system goes to the ground state as expected. In fact, we found that 
at zero temperature the mean valence is determined by the parameter vq and we approximated 
the first order correction, showing the dependence on the coupling constants of the model. While 
these results were expected on general grounds, the mapping used here simplified the problem and 
allowed a quantitative analysis. 

What emerged from the study of the average distribution for the valence is that, if the vertex 
valence term dominates (gy 3> g^), the loop term corrections to the average valence of the ground 
state are suppressed at T = 0. We found the dependence on the coupling constants explicitly. 
This result is confirmed by the study of the susceptibility. Thanks to the mean field theory 
approximation, we found the contribution of all loops to the susceptibility and showed that the 
susceptibility function tends to the Ising one when T approaches zero. We then applied this 
procedure to the degenerate case, and showed that the degeneracy does not change the average 
valence at low temperature but only the speed with which this ground state is reached. As a final 
remark, we stress that the vertex valence is an important quantity in the model. In fact, as shown 
in |19j using the Lieb-Robinson bound, the speed with which information can propagate on graphs 
is bounded by a valence-dependent quantity. For this reason, as the temperature drops, the speed 
of the emergent light field must drop with the valence. 

We would like to stress that we assumed that the system was at equilibrium with an external 
bath. The problem of the required external bath in the model has been studied in [16], where 
additional degrees of freedom (bosonic particles on vertices) were introduced to let the system 
thermalize and reach an equilibrium distribution. 
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